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Abstract. For the third g-Bessel function (first introduced by F. H. Jackson, later 
rediscovered by W. Hahn in a special case and by H. Exton) we derive Hansen- 
Lommel type orthogonality relations, which, by a symmetry, turn out to be equivalent 
to orthogonality relations which are (/-analogues of the Hankel integral transform pair. 
These results are implicit, in the context of quantum groups, in a paper by Vaksman 
and Korogodskii. As a specialization we get (/-cosines and (/-sines which admit q- 
analogues of the Fourier-cosine and Fourier-sine transforms. We get also a formula 
which is both an analogue of Graf's addition formula and of the Weber-Schafheitlin 
discontinuous integral. This is a corrected version of a paper which appeared in 
Trans. Amer. Math. Soc. 333 (1992), 445-461. 



1. Introduction 
Several possible (/-analogues of the Bessel function 

(11) J (x) ■= f ^ (X/2)a+2k 
1 > a[ '■ £j fc!r(a + fc+l) 

have been considered in the literature. Best known are two related g-Bessel func- 
tions denoted Ja\x; q) and J&\x; q) by Ismail [10], but first introduced by Jackson 
in a series of papers during the years 1903-1905 (see the references in [10]) and also 
studied by Hahn [7]. A third g-Bessel function was introduced by Hahn [8] (in a 
special case; we thank G. Gasper for this reference) and by Exton [3], [4, (5.3.1.11)] 
(in full). In Exton's notation: 

(12) CJq- x) - (1 ~ ^ T qHk ~ 1)/2 {qa+k+1; q) °° { ~ X(1 - ^ 
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This paper appeared in Trans. Amer. Math. Soc. 333 (1992), 445-461. In the present version 
the statement and proof of Proposition 4.2 and the proof of Proposition A.l have been corrected. 
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where the (/-shifted factorials are defined by 



k-1 



(1.3) 



(a;q)k--=Y[(l-aq J ), fc = l,2,...; (a;g) :=l; 

3=0 



(a;q) 



oo 



lim (a;g) fc , 



q\ < 1. 



Hahn [8] considered the case a = of (1.2). They obtained these functions as 
the solutions of a special basic Sturm-Liouville equation, by which they could also 
derive the following (/-analogue of the Fourier-Bessel orthogonality relations: 



where fj,i, /i2, • • • are the roots of the equation C a (q; fx) = 0, and where the (/-integral 
is defined by 



By specialization to a = ±1/2, Exton obtained similar orthogonalities for q- 
analogues of the sines and cosines. So some of the harmonic analysis involving 
Bessel functions, sines and cosines has been extended to the g-case. However, q- 
analogues of the Fourier-cosine, Fourier-sine and Hankel transforms were missing 
until now (except for a (/-Laplace transform with inversion formula, cf. Hahn [7, §9] 
and Feinsilver [5]). 

Exton [4, (5.3.3.1)] also generalized the generating function 




i^j, a > -1, 




oo 



(1.4) 




k=- 



oo 



(cf. Watson [15, §2.1 (1)]) to the case of his (/-Bessel functions (1.2): 



Ov: 



(1.5) 



e q ((l-q)t) E q (-(1 -q)t~ 1 x)= £ t n C n (q;x), 



n= — oo 



where 



(1.6) 




z\ < 1, 



and 



(1.7) 
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are (/-analogues of the exponential function (cf. Gasper and Rahman [6, (1.3.15), 
(1.3.16)]). 

Recently, Vaksman and Korogodskh [14] gave an interpretation of the g-Bessel 
functions (1.2) as matrix elements of irreducible representations of the quantum 
group of plane motions. Their paper, which does not contain proofs, implicitly 
contains some new orthogonality relations for the functions (1.2). In particular, 
the unitariness of the representations implies a (/-analogue of the Hansen-Lommel 
orthogonality relations 

oo 

(1.8) S nm = ^ J k+n(x) J k+m (x), n,meZ 

k= — oo 

(cf. [15, §2.5 (3), (4)]). Furthermore, the Schur type orthogonality relations for 
matrix elements of irreducible unitary representations which are square integrable 
with respect to a suitable Haar functional, yield a (/-analogue of Hankel's Fourier- 
Bessel integral 

/•oo />oo 

(1.9) f(x)= J a (xt)( J a (ty)f(y)ydy)tdt, 

Jo Jo 

cf. [15, §14.3]. 

It is the purpose of the present note to state these two types of orthogonality 
relations for the functions (1.2) explicitly, to show that the first type is immediately 
implied by the generating function (1.5), and to rewrite the first type as the second 
type by use of a simple, but possibly new, symmetry for the functions (1.2). This 
will be done in §2. In §3 we will show that the second type of orthogonality is, on the 
one hand, a limit case of the orthogonality for the little (/-Jacobi polynomials, and, 
on the other hand, allows the Hankel transform inversion formula as a limit case 
for q t 1. In §4 we will generalize the two orthogonality relations to two equivalent 
formulas, which are respectively the g-analogues of Graf's addition formula and the 
Weber- Schafheitlin discontinuous integral. The special cases of the (/-Fourier-cosine 
and sine transforms will be the topic of §5. No material from §4 is needed in this 
section. Finally, Appendix A will contain rigorous proofs of two limit results. 

In a recent paper by Rahman [13], where a (/-analogue of the Fourier-Bessel 
orthogonality for Jackson's (/-Bessel functions is discussed, the author states in his 
concluding remarks that Jackson's g-Bessel functions have probably nicer properties 
than those of Exton. However, the results obtained in [14] and in the present 
paper might suggest that the Hahn-Exton functions are more suitable for harmonic 
analysis, both within and without the context of quantum groups. Future research 
will help to clarify the merits of the various types of g-Bessel functions. 

In this paper we will not preserve Exton's notation C a (q;x) in (1.2), but state 
the results in terms of the g-hypergeometric function 

(1-10) 101 (0; w- q, z) := £ 1 / \ . 

{w,q)k(q]q)k 
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Our motivation is that (i) the r cj) a notation is fairly well known nowadays and clari- 
fies the position of these (/-Bessel functions among other g-hypergeometric functions; 
(ii) g-analysis should not depend too much on the q = 1 case, so scaling factors 
making the limit transition q |" 1 easy should not be hidden in the definitions of 
(/-special functions; and (iii) the classical definition (1.1) of Bessel functions is very 
natural in the context of the generating function (1.4) and for analysis concentrat- 
ing on J , but it is less fortunate when the focus is on another Bessel function of 
fixed order (cf. J_i/ 2 (x) versus cosx), so one should be very careful before fixing 
the definition and notation of a (/-Bessel function. Therefore, the notations for q- 
Bessel functions, (/-cosines and (/-sines in §§3 and 5 should be considered as ad hoc 
notations, which are only used locally in this paper for clarifying the analogy of the 
functions with the q = 1 case. 

While this Introduction until here describes what was known to us in 1992, it 
has later been observed by M. E. H. Ismail (see for instance p. 184 in [16]) that this 
third type of g-Bessel function was already observed by Jackson on p. 201 of his 
paper [17]. His notation and definition is 

T m— -\\r n r(r-l) ( (1 ~ ff) A) n+2r 

J[.l(A).-L,(-l)p' (p2 ;f)2) „ +r(p2;p2)/ 

Therefore, what we called the Hahn-Exton g-Bessel function is now often called 
Jackson's third q-Bessel function. 

2. Symmetry and orthogonality for q-Bessel functions 

We will always assume that < q < 1. The general (/-hypergeo metric series is 
defined by 

(a 1 ,...,a r] q) k ((-l)U« k - 1)/2 Y~ r+1 z k 
(b 1 ,...,b s ;q) k (q;q) k 

where the (/-shifted factorial is defined by (1.3) and 

(ai, . . . , a r ; q) k := (a x ; q) k (a 2 ; q) k • • • (a r ; q)k- 

The upper and lower parameters in the left-hand side of (2.1) may also be written 
on one line as ai, . . . , a r ; bi, . . . , b s . The power series, in the non-terminating case 
of (2.1), has radius of convergence oo, 1 or according to whether r — s < 1, = 1 
or > 1, respectively (see [6, Ch.l] for further details). Thus the defining formula 
(1.2) for the Hahn-Exton g-Bessel function can be rewritten as 

C a (q; x) := (1 ~ 9) " (9 " +1;9)o ° 10! (0; (/"+ 1 ; q, x(l - qf) 



(2.1) r 4>s 



ai, . . . , a r 
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and the (/-exponential functions (1.6), (1.7) can be written as i0o(O; — ; q, z) and 
o0o(-; -; Q, ~z), respectively. 

Our object will be the Hahn-Exton g-Bessel function written as the g-hypergeometric 
series (1.10). It is well-defined for z, w G C with w outside {1, g -1 , q~ 2 , . . . }. These 
singularities can be removed by multiplication by (w;q)oo: 



(2.2) (w;q) ool (j> 1 (0;w;q,z) = ^2 



(-l) k q k( - h 1 ^ 2 (q k w;q) QO z h 



k=0 

Proposition 2.1. The series in (2.2) defines an entire analytic function in z,w, 
which is also symmetric in z, w: 

(2.3) {w,q)oo i(f>i{0;w;q,z) = (z^q)^ i0i(O; z; q,w). 
Both sides can be majorized by 

(2.4) (-M;?)oo(-H;9)oo- 

Proof. Substitute for (q k w;q)oo in (2.2) the o</>o series given by (1.7): 

(_!)* g i('-i)/2 w l q k(k-i)/2 z k 



(w;q)oo i0i(O;«;;g,z) = $^9*' 



(?;?) 



k=0 1=0 

The summand of the double series can be majorized by 

q l(l-l)/2 j^ji q k(k-l)/2 | z jfc 



(?;?)fc 



(9; 9) 



(9;s)fc 



Thus the double sum converges absolutely, uniformly for z, u; in compacta, and it 
is symmetric in z and w. □ 

Remark 2.2. Formula (2.3) is a limit case of Heine's transformation formula 



(2.5) 



201 



a, b 



(6,az;g) c 
(c,z;g)o< 



2</>l 



c/6, ^ 



5 9,6 



(cf. [6, (1.4.1)]). Indeed, first let b — > 0, then replace 2 by z/a and let a — > oo. At 
least formally, by termwise limits, we obtain (2.3). 

Remark 2.3. For w := q x ~ n (n = 1, 2, . . . ) we interpret the left-hand side of (2.2) 
by the series at its right-hand side. Then the first n terms vanish, so the summation 
starts with k = n. When we make the change of summation variable k = n + /, we 
obtain 
(2.6) 

(q 1 - 71 ; q)oo i^i (0; q X ~ n ] q, z) = (-l) n ^" 1)/2 z n (q n+1 ; q)^ ^(0; q n+1 ; q, q n z) 

for n E Z. (The case n < follows from the case n > of (2.6) by changing z into 
q~ n z.) 
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Remark 2.4. Because of (2.6), the behaviour of the two equal sides of (2.3) as 
\w\ — ?■ oo (cf. (2.4)) drastically improves when w runs over the values (/ 1_n , n = 
1,2, For such w we can majorize these expressions by 

<z n(n - 1)/2 M n (-M;<z)oc(-^)oo. 



We will now restate Exton's generating function (1.5) in terms of the notation 
(2.2), and also give the short proof, for reasons of completeness. 

Proposition 2.5. For z,i G C such that < \t\ < there is the absolutely 

convergent expansion 

(2-7) 



(tz;q)oc „fr^ ' (q; q) 



oo 



(2.8) = t n z ni -f^ 1 M0;z 2 ;i,Q n+1 )- 



Proof. Expansion of the left-hand side of (2.7) gives 

~ (-l) fc gfc(fc-D/2 t l-k z l+k _ ~ ~ (g i+l. g)oo( _ 1)fcg fc(fc-l)/2^-fc^+fc 

which is an absolutely convergent double sum for z, t G C such that < \t\ < | ^ | 1 . 
Now pass to new summation variables k, n by substituting I = k + n. This yields, 
by substitution of (2.2), the right-hand side of (2.7). Formula (2.8) follows by the 
symmetry (2.3). □ 

In §3 we will show that (2.8) is a (/-analogue of an infinite integral of Weber and 
Sonine. Replace in (2.7) t by t~ x and multiply the new identity with the original 
identity. The resulting formula 

oo oo / n+1 \ 

( 2 - 9 ) 1= E E f ■ \ q)o ° iM0;q n+1 ;q,z 2 ) 

\Qi Q)oo 



n= — oo m= — oo 



00 ± (r\ „m+l „ „2\ 

xz — - i0i O g ;q, z ) 



I oo 



is absolutely convergent for z, t G C such that t ^ and \z\ < \t\ < So 
equality of coefficients of equal powers of t at both sides yields a (/-analogue of the 
orthogonality relations (1.8): 
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Proposition 2.6. For \z\ < 1 and n,m G Z we have 

00 f„n+fc+l.„\ 
In 1n \ ^fc+n W > V/oo i In. „n+fc+l. _ „2\ 

(2.1UJ 2^ — r - 1 i0h u ; ( ? ;<?^ J 



(?;?) 

and 



/m+fc+l. „\ 

7 v i0i(O;g T T j = o. 



n rn 



:3C 



(2.11) £ z k+n{ p4^iM0;^,Q,Q n+k+1 ) 

k=-oo ^ q) °° 

^ k+m ( Z 'iQ)co , /„ 2 „ m+fc+l\ r 

xz T r — i0i (0; z ;q,q T T ) = o nm , 

where the sums at the left-hand sides are absolutely convergent, uniformly on com- 
pact subsets of the open unit disk. 

Formula (2.11) follows from (2.10) by the symmetry (2.3). In §3 we will show 
that (2.11) is a g- version of HankePs Fourier-Bessel integral (1.9). 

Remark 2.7. Analogous to Proposition 2.5 there are the two generating functions 
for the g-Bessel functions J&\x;q) and J&\x;q) of Jackson and Ismail [10]: 

(2.12) e q (tz)e q (-t- 1 z) = 



(tz; g)oo (-t~ 1 z;q) c 
j n+1 ;q 



OO / n -|-l_ \ 

E,n n KQ 'iQjco i / n „ n+1. „ 2\ I i ^ |.| ^ I „|-1 
t z ^ 2 0i(O,O;g ;q,-z ), \z\ < \t\ < \z\ , 



loo 

li — — {JU 

and 

(2.13) E q (tz) E q (-t~ l z) = (-tz; g)^ (t^z; g)^ 

00 „n(n-l)/2 ( n n+l. „\ 

In a way similar to (2.10) we can now derive the biorthogonality relations 

00 ( n n+k+l. n \ 

(2.14) 2^ r S 2<Pi(0,0;g T T ; g, -z ) 

/„m+fc+l. „\ 

v r fc+m _(m+fc)(m+fc-l)/2 W ; ( i)oo il . _m+fc+l._ n m+k 2\ _ r 

x ^ g — 7 — s o<m — >g ,q,—q z) — o nm , 

{q;q)oo 

valid for |z| < 1 and n, m G Z. The case n = m of this result goes back to Jackson 
(see also [7, §8]). Formula (2.14) can be rewritten in several ways by substitution 
of the transformations 

o^i (-; c; g, cz) = (z; g)^ 2 0i(O, 0; c; g, z) = — i0i(z; 0; g, c). 

(c;g)oo 

However, this will not transform (2.14) into an orthogonality; it remains a biorthog- 
onality. 
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3. Some limit transitions 
Jacobi polynomials tend to Bessel functions: 

P:::; 1 ( :i - x 2 /(2N 2 J)_ = 2i?i( _ n ^ nN + a + p + 1 . a + 1 . x 2 /(4Ar2)) 



"> Fi(-; a + 1; -(Ax/2) 2 ) = (\x/2)~ a T(a + 1) J a (Xx), 



iV-s-oo 



where un/N tends to A for N — > oo. When this limit transition is applied to the 
formula which recovers a function from its Fourier- Jacobi coefficients, we obtain, 
at least formally, Hankel's Fourier-Bessel integral (1.9). 

The (/-analogue of this limit transition starts with the little q-Jacobi polynomials 

p n (x;a,b;q) := 2 <M<?~ n , abg n+1 ; ag; q, qx), 

which satisfy orthogonality relations 



loo 
) oo 



Sr. 



(31) (qa,qb;q)oo yv )( a k - a b- a) (aa) k iS^lllk 
(3 - 1} {q,q*ab,qU^ PnPm)[q ' a ^ q)[qa) (q^b;q\ 

= (ga) n (1 - gab) (gb,g; g) n 
(1 - q 2n+1 ab) (qa,qab; q) n 

where < a < (/ _1 , b < q~ l (see Andrews and Askey [2]). 
It is clear that 

p N - n (q N x; a, b; q) = 2 MQ~ N+n , abg N ~ n+1 ; aq; q, q N+1 x) 

tends formally (termwise) to i0i(O; aq; q, q n+1 x) as N — > oo. (See Proposition A.l 
for a rigorous proof of this limit result.) Also, when we replace in (3.1) n, m, k by 
N — n, N — to, N + k, respectively (so the sum runs from —N to oo), and when 
we let N — > oo, we obtain as a formal (termwise) limit the orthogonality relations 
(2.H). 

In order to see that (2.11) is a (/-analogue of Hankel's Fourier-Bessel integral 
(1.9), rewrite (2.11) as the transform pair 



(3.2) { 



oo / 2a+2 2\ 

k^oo (Moo 



oc 



f(q k ) = E ? (*+»)(a+D ^ a 2 +2; 2 f )o ° ^ 2a+2 ; <?> <z 2fc+2n+2 ) g(g n ), 



where /, (7 are L 2 -functions on the set {g k \ k £ Z} with respect to counting measure. 
With the ad hoc notation 

a ( n 2a+2. 2\ 

/q n\ j /„. J2\ . * \H > g )oo 1 / n 2a+2. "2 „2 J2\ 

(6.6) J a (z;g):= 2 _ 2 i0i(U;g ^ ;g ,g z ) 

(g ? </ joo 
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(different from Rahman's proposal in [13, (1.13)]) and with the replacement of 
/(<?*% g(q n ) by q k f(q k ), q n g(q n ) this becomes 



(3.4) 



Now observe that 



9(q n )= E q 2k Uq k+n ;q 2 )f(q% 

k= — oo 

oo 

f(q k )= E q 2n Jc(q k+n ;q 2 )g(q n )- 



J Q ((1 - q)z; , 2 ) = ( ^ + ; } "7" 101(0; , 2 "+ 2 ; , 2 , (1 - , 2 ) 2 (qz/(l + q)f) 
T q 2(a + 1) 

converges, for q f 1, to the Bessel function 

= F ( a + 1 ) o^i (-; a + 1; -^ 2 /4), 

where we used Proposition A. 2 and the fact that 

(1 - <Z 2 ) a (<7 2Q+ W)oo 1 . 1 



(q 2 -^ 2 )^ r g2 (« + i) r(a + i) 

as 5 t 1, cf. Andrews [1, Appendix A], Koornwinder [11, Appendix B]. We can 
apply this to (3.4) when we let q ~\ 1 under the side condition that lo ^^~ <? ' > e 2Z. 

For such q we can replace q k , q n in (3.4) by (1 — q) 1 ^ 2 q k , (1 — q) l / 2 q n , and next 
/((l — qy/ 2 q k ), g((l — g) 1 / 2 q ,n ) by /(g fe ), g{q n )- With the (/-integral notation 



/ h(z)d q z:=(l-q) V %>V 



j = -oo 

(3.4) then takes the form 



/•OO 

gW = f(x)J ol ((l-q)\x;q 2 )xd q x, 
Jo 

/•oo 

f(x)= / s(A) J a ((l-q)Xx;q 2 )Xd q X, 
Jo 



where A in the first identity and x in the second identity take the values q n , n G Z. 
For q "f 1 we obtain, at least formally, the Hankel transform pair 



/•OO 

fi'(^) = / f(x) J a (Xx)xdx, 
Jo 

I* OO 

/(x)= / y(A) J a (Ax)ArfA, 
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which is equivalent to (1.9). 

In order to find the classical formula corresponding to (2.8), we rewrite (2.8) first 
in terms of the notation (3.3): 

( n a+l-t 2\ °° 
W > Q } °° n=-oo 

For log(l — q)/\ogq £ Z this can be rewritten as 



(l + g)*I>((a+l + *)/2) 



I>((a + l-*)/2) 

Formally, as g f 1) this yields 

2 t r((q + l + t)/2) r°° 
T((a + l-t)/2) 



/•OO 

/ x t J Q ((1 - q)x;q 2 ) d q x. 
Jo 



which formula, valid for —1/2 > 9fc£ > — 3?ct — 1, goes back to Weber and Sonine 
(cf. [15, 13.24 (1)]). 

4. A (/-ANALOGUE OF GRAF'S ADDITION FORMULA 

In this section we will generalize the considerations which led to the orthogonality 
relations in Proposition 2.6. The resulting formula will turn out to be a (/-analogue 
of Graf's addition formula and, at the same time, of the discontinuous integral 
of Weber and Schafheitlin. In a final remark we will point out that the Graf 
type addition formula part can also be done for the (/-Bessel functions of Jackson 
and Ismail, by similar methods, which go back to Heine, the founding father of 
(/-hypergeo metric series. 

A formula more general than (2.9) can be derived by expanding the expression 

^ 4 -g (xs- 1 t;q) 00 (yt- 1 ;q) 00 



(yt;q) QO (xst 1 ;q) 00 

as a Laurent series in t (\sx\ < \t\ < |y| _1 ) in two different ways. On the one hand, 
(4.1) can be expanded by twofold substitution of the (/-binomial formula (cf. [6, 
(1.3.2)]) as: 

1 (po(s~ 1 xy~ 1 ] -; q, yt) 1 (p (s~ 1 yx~ 1 ; -; q, xst' 1 ) 

= y y (s- 1 yx- 1 ;q) k (s- 1 xy- 1 ;q) l gk ^ f 

= f t n y n y (s-yx-^q), (s-ixy-^q) n+k k ^ 

n^oo fc=ma tf- n ,0} («^)* 
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where we substituted / = k+n. Since the inner sum in the last part can be expressed 
in terms of a 291 series (depending on the sign of n), we obtain for \sx\ < \t\ < 
the identity 

(4-2) 

{xs~H, yt' 1 ; q)oo ^ xn _. n (s^xy' 1 ; q) r 



(yt,xst 1 ;q) c 



E*v 



+ E * 



n=0 



n s~ n x 



291 



s x yx 1 ,q n s 1 xy 1 



7 n+l 



; q, sxy 



291 



s 1 xy 1 ,q n s x yx 

n -n+l 



; q, sxy 



(«;«)-» 

We might have written (4.2) as a sum over n running from — oo to oo with summand 
given for all integer n by the left-hand side of (4.3) below, where, for n < 0, an 
interpretation of the 291 similar to our convention in Remark 2.3 is used. Thus the 
analogue of (2.6) becomes: 



(4.3) t n y' 



(s 1 xy 1 ,g 1+n ;g) c 
(q n s- 1 xy- 1 : q;q) 0i 

t n s- n x- n 



291 



s x yx 1 ,q n s x xy ' 



,1+n 



; q, sxy 



(s 1 yx 1 ,q 1 n ; 



291 



s l xy 1 ,q n s x yx J 

n l-n 



; 3, sxy 



(q n s 1 yx 1 ,q;q) c 

Here we should exclude the case that both 1 + n and n + q log(s~ 1 xy~ 1 ) are non- 
positive integers, or that both 1 — n and — n + q log(s~ 1 yx~ 1 ) are non-positive 
integers. 

Note that (4.2) reduces to (2.7) in the special case x = (and also, in view of 
(2.6), for y = 0). Formula (4.2) is a (/-analogue of (1.4) with z and t replaced by 
2(y — s~ 1 x) 1 / 2 (y — sx) 1 ^ 2 and t(y — s~ 1 x) 1 ^ 2 (y — sx) -1 / 2 , respectively. 

On the other hand we expand (4.1) by twofold substitution of (2.7): 



(4.4) 



(as t;q) 00 (yt 1 ;q) 

(yt;q)oo (xst-^q) 



= E E ^ 

n= — oo k= — oo 



X y 



„+k (q n+k+1 ; g)oo . n+fc+1 _ 2 \ ^k (Q fc+1 ; g)c 



i0i(O; 9 n+ * +1 ;3,ir)*' 



1 1 (O;g fc+1 ;g,x :i ), 



(9; 9)00 (g;g)< 
which generalizes (2.9). When we compare coefficients of equal powers of t in (4.2) 
and (4.4), we obtain: 

Proposition 4.1. For \sxy\ < 1 we have 
(4-5) 

£ (gW+fc+ ; g) °° i0i (0; «" +fc+1 ; «, 2/ 2 ) ^ ^ (0; ? , x 2 ) 

k=-oo KQ,Q)oo [q,Q)c 



V 



291 



s- n x~ n 



(q;q) n 

(s~ 1 yx~ 1 ;q) 



s x yx 1 ,q n s x xy 1 

n n+l 



q, sxy 



if n > 0, 



(?;?)- 



2</>l 



s 1 xy ,q n s 1 yx 1 

n-n+l 



q, sxy 



if n < 0. 



12 



TOM H. KOORNWINDER AND RENE F. SWARTTOUW 



Formula (4.5) is a (/-analogue of the addition formula 

n/2 



(4.6) J2 s k Jn+M Mx) 



k= — oo 



y — s x 
y — sx 



Jn[v(y- s lx )(y ~ sx ) 



due to Graf (cf. [15, §11.3 (1)]). The special case n = of (4.5) is a (/-analogue 
of Neumann's addition formula for Bessel functions Jo (cf. [15, §11.2 (1)]). In the 
special case x = y, s = 1 the right-hand side of (4.5) becomes y n 5 n fi, so then (4.5) 
reduces to the orthogonality relations (2.10). 

When we apply the symmetries (2.5) and (4.3) to the left-hand side of (4.5) and 
also apply (2.3) to the right-hand side of (4.5), and replace n by n — m and next k 
by k + m then we obtain an equivalent identity: 

Proposition 4.2. For \sxy\ < 1 and x 2 , y 2 ^ q~ l (I = 0, 1, 2, . . . ), we have 
(4-7) 

g k y n + k (^ 2 '^)°° A (C\.„2.„ „n+k+l\ „m+k O^i <Z)c 



0i(O;y ;<?,<? 



)x r 



(<?;<?)c 



101 (0; x ;q,q 



m+fc+l 



s~ m y n ~ m 



(g~ 1 a;y~ 1 ,y 2 ;g) c 
(sxy.q-.q)^ 

(g~ 1 yx~ 1 ,x 2 ;g) c 
(sxy.q-.q)^ 



291 



291 



qsx x y,sxy n _ m _ x _ x 

2 ;q,q s xy 

y 

qsxy^.sxy m _ n _ x _{ 
2 ;q,q s yx 
x z 



where the first expression on the right-hand side can be chosen if s x xy 1 and 



f l ~ m s~ 1 xy~ 1 are not equal to a non-positive integer power of q, and where the sec- 
ond expression on the right-hand side can be chosen if s~ 1 yx~ 1 and q m ~ n s~ 1 yx~ 1 
are not equal to a non-positive integer power of q. Thus in the generic case both 
expressions on the right-hand side are equal. 

See Appendix B for the detailed proof. Observe the special case x = y of (4.7): 



(4.S 



k= — oo 



fc ^n+k ( z2 'i q)oo i (n.. J2. „ „n+k+l\ ^m+k Q)oo a / n . Ji. r , r ,m+k+l\ 

s z ; — i0i (0; z ;q,q T T )z — i(pi{U; z ;q,q T T ) 



(q;q) 



oc 



(q;q) 



■oc 



s -m z n-m 



s -n z m-n 



(s-\z 2 ;q) 
(sz 2 ,q; q) c 
(s-\z 2 ;q), 
(sz 2 ,q;q) a 



291 



291 



(/S, SZ n—m — 1 

Z A 

qS 'f;q,q m - n s-' 
z 2 



where \sz 2 \ < 1, z 2 is not a non-positive integer power of q and the expression on 
the right-hand side has to be chosen similarly as in Proposition 4.2. This can be 
considered as a kind of Poisson kernel for the orthogonal system with orthogonality 
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relations (2.11). Inspection of the right-hand side of (4.8) shows that this kernel is 
positive if < z < 1 and 1 < s < min{(/ _1 , z~ 2 }. 

Let us look for the classical analogue of formula (4.7). In (4.7) first replace q by 
q 2 , then x,y,s by q a+1 , 9 /3+1 , <Z _7_1 , respectively. Then, with the notation (3.3), 
formula (4.7) can be rewritten as 



£ q k ^ur +k ;d 2 )M<i n+k ;<i' 



k= — oo 



n/3 m( 7 -/3-l) W )c 



x 2 0l(/- Q - 7+1 , q a +^ +1 ; q 2 ^ 2 - q 2 , ? 2n-2m+a-/3+ 7 +l ); 



/„/3-a+7+l n 2a+2. n 2\ 
ma n(y-a-l) V± ij± i H )ce 

q q (q a+ ^+\q 2 ;q 2 )oo 



x 2<M<Z Q ~ /3 ~ 7+ \ q a+f) -i +1 ; q 2a+2 ; q 2 , q ^-2n+p-a+ 1+ i^ 



where 3? (a + /3 — 7 + 1) > 0, a, /3 are not equal to a negative integer, and where 
we can choose the first expression on the right-hand side if (a — f3 + 7 + l)/2 and 
n — m + (a — (3+ / y + l)/2 are not a non-positive integer and we can choose the second 
expression on the right-hand side if (/3 — a + 7 + l)/2 and m — n+(f3 — a + 7 + 1)/2 
are not a non-positive integer. 

Now replace <j» fc by q k (l — 5) (with ^f^y^ £ Z) and let m,n depend on g such 
that, as q | 1, q m tends to a and q n tends to 6. Depending on whether b < a 
or 6 > a, make the formal limit transition q "f 1 in the first or second identity, 
respectively. Then, for 3?(a + /3 — 7 + 1) > 0, 3?7 > —1, we obtain the discontinuous 
integral of Weber and Schafheitlin: 

poo 

2 7 / x -1 J a {ax) Jf3(bx) dx 
Jo 

' Q 7-/3-i r((q + ft - 7 + l)/2) 
r((a-/3 + 7 + l)/2)r(/3 + l) 

x 2 F 1 ((/3-a-7 + l)/2,(a + /3- 7 + l)/2;/3 + l;6 2 /a 2 ) if 6 < a, 
a Q r((a + g - 7 + l)/2) 

r((/3-a + 7 + l)/2)r(a + l) 

x 2 F 1 ((a-/3-7 + l)/2,(a + /3-7 + l)/2;a + l;a 2 /^ 2 ) if a < 6, 



(cf. [15, §13.4 (2)]). Note that the two analytic expressions at the right-hand side 
are no longer equal, as they were in the g-case. 
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Remark 4.3. Analogous to (4.2) we have 
(xs^t^q)^ (-xst- 1 -^)^ 



(4.9) 



(yt;q)oo(-yt 1 ;q) 00 

X 201 



V- f n n ( S J?® 1 i g" +1 ;g) 
2-^ ^ (nns-lrro.-l 



oo 
oc 



q n s 1 xy 1 ,sxy 1 2 

qn +i '^-y 



(q n s x xy 1 ,g;g) c 

, \y\ < \t\ < M" 1 - 



The case s = 1 of (4.9) goes back to Heine [9, p. 121] (see also Hahn [7, §8]). Its 
special cases (x,y,s) = (0,2,1) and (—2,0,1) are the formulas (2.12) and (2.13). 
Like (4.2), formula (4.9) is a (/-analogue of (1.4) with z and t replaced by 2(y — 
s~ 1 x) 1 ' 2 (y — sx) 1 / 2 and t(y — s~ 1 x) 1 ^ 2 (y — sx) -1 / 2 , respectively. Similarly to (4.5) 
we obtain from (2.12), (2.13) and (4.9) that, for \y\ < 1: 



(4.10) y' 



(s x xy 1 ,g n+1 ;g) 00 
(q n s- 1 xy- 1 ,q;q) (> 



29l 



^r^^ 00 201 (0, 0; q k+n+1 ; q, -y 2 ) x k q k ^' 2 ^ + \^°° o l( - ; q , - q * x 2 



q n s x xy x ,sxy 1 2 

qn+ i '^-v 

fc+1. 



= E 

k= — oo 



,k+n 



X 



(g;g) 



The special case n = of this formula is the limit case v \. of Rahman's addition 
formula [12, (1.10)]. Like (4.5), formula (4.10) is a g-analogue of Graf's addition 
formula (4.6). The special case x = y, s = 1 of (4.10) gives the biorthogonality 
relations (2.14). 

5. q- Analogues of the Fourier-cosine and Fourier-sine transform 
Put 



(5.1) 
(5.2) 
and 

(5.3) 
(5.4) 



cos(2;g 2 ) :=i0i(O;g;g 2 ,g 2 2 2 ) 

(-l) fc q k ( k+1 ) z 2k 



=E 



k=0 

(g 2 ;g 2 )c 
(q;q 2 )o 



(g;g) 



2 k 



z^ 2 J_ 1/2 (z-q 2 ) 



sin(2; q 2 ) :=(1 - g)" 1 z ^(0; g 3 ; g 2 , g 2 z 2 ) 
(-l) k q k ( k+1 1 z 2k+1 



k=0 
,2. „2 



(g; g) 



2fc+l 



(g ; g )oo 1/2 t / 2\ 
(g;g 2 )oo 
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Here we have used the notation (3.3). (The functions introduced above should 
not be confused with the functions cos 9 and sm q considered in [6, Exercise 1.14].) 
Clearly we have the formal (termwise) limits 



cos((l — q)z; q 2 ) — Ycosz and sin((l — q)z; 



sin ^ 



as q t 1. By Proposition A. 2 these limit transitions hold pointwise, uniformly on 
compacta. When we substitute (5.2) or (5.4) in (3.4) and replace f(q k ), g(q n ) by 
q~ k ^ 2 f(q k ), q~ n ^ 2 d(q n ), we obtain the transform pairs 



(5.5) 



9(q n ) 



f(q k ) 



{q;q 



2\ 00 



(q;q 2 



{q 2 ;q 2 )c 



cos(q k+n ; q 2 ) 

] or ) f(q k ), 

sm(q k+n ; q 2 ) 

cos(q k+n ; q 2 ) ' 

or \ g{q n ). 

y sin(g fc+n ; q 2 ) ^ 



E o l 

fc= — oo 
oo 

E 



The transformations / i-> g and g i— y f of (5.5) establish an isometry of Hilbert 
spaces: 

oo oo 

E ? fc i/(^)i 2 = E ? n i^ n )i 2 - 

fc= — oo n= — oo 

Now let q p 1 under the side condition that lo ^^~ 9 - > e 2Z. Replace g n in (5.5) 

by (1 - g) 1 / 2 ^, (1 - g) 1 / 2 ^ and then /((l - g) 1/2 g fc ), #((1 - g) 1/2 g n ) by /(g fc ), 
g(q n ). Then (5.5) takes the form 



<7(A) 



/(*) 



(1 + g) 1 / 2 

r>(i/2) 



poo 

/ /( 

^0 



x) 



cos((l — q)Xx; q 2 ) 



or 



(1 + 9) 



1/2 /-oo 



i> a/2) y 



<?(A) 



sin((l — q)Xx; q 2 ) 
cos((l — q)Xx; q 2 ) 



or 



sin((l — q)Xx; q 2 ) 



dqX, 



d q X. 



Formally, as q 1 1, we obtain the classical Fourier pairs 



and 



/•oo poo 

g{X) = a/2//t / /(x) cos(Ax) dx, f(x) = a/2/ 71 " / g(X) cos(Ax) dX 
Jo Jo 

pOO pOO 

g(X) = / f(x) sin(Xx) dx, /(x) = ^2pH \ g(X) sin(Ax) dX. 

Jo Jo 



With the notation 



(D q f)(z) = D q , z f(z) := 



/(*) ~ f(qz) 
(l-q)z 
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for the (/-derivative, we obtain from (5.1) and (5.3) that 

(1 - q) D QjZ cos(z; q 2 ) = -q sm(qz; q 2 ), 
(1 - q) D qjZ sm(z; q 2 ) = cos(z; q 2 ). 



Hence 



(±-q) 2 (D 2 q f)(q- 1 z) = l [ 



-q 2 X 2 f(z) iff(z) = cos(Xz;q 2 ), 
-q\ 2 f(z) iff(z) = sm(Xz;q 2 ). 



So the two systems of functions z h-> cos(q n z) (n G Z) and z i— >■ sm(q n z) (n G Z) 
have disjoint eigenvalues with respect to the operator which sends / to the function 

^(1- Q) 2 {Dlf){q^z) =qz- 2 {qf(q- 1 z) - (1 + q) f(z) + f(qz)). 

This operator has also the selfadjointness property 

oo oo 
A:= — oo A:= — oo 

for /, g of finite support on {q k \ k G Z}. 

Observe that the (/-deformation of d 2 /dx 2 considered above yields a symmetry 
breaking. The two-dimensional eigenspaces of d 2 /dx 2 are broken apart into one- 
dimensional eigenspaces. Therefore it does not seem to be very useful to consider a 
g-exponential built from the functions defined by (5.1) and (5.3). Any linear com- 
bination f(z) of cos(\z;q 2 ) and sin(Xz;q 2 ) will no longer satisfy an eigenfunction 
equation 

(5-6) (l-q) 2 (D 2 q f)(q- 1 z) = ^f(z), 

while the nice function 

f(z) := ^(O; -q 1 / 2 - q 1 / 2 , ±iq 3/4 \z) = cos(Xz; q 2 ) T V /4 smtf^Xz; q 2 ) 

(cf. Exton [4, 5.2.2.1]), which satisfies (5.6) with \i = —q 2 X 2 : no longer remains 
within the spectral decompositions implied by (5.5). 

Appendix A. Rigorous proofs of some limit results 

In this Appendix we will give proofs of the limit transitions from little (/-Jacobi 
polynomials to (/-Bessel functions and from (/-Bessel functions to ordinary Bessel 
functions. 
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Proposition A.l. For < a < q~ l , < b < q~ l we have 

lim 2 4>i(q~ n , Q n+1 ab; qa; q, q n x) = i0i(O; qa; q, x), 

n— S-oo 

uniformly for x in compact subsets of C. 
Proof. Put 

R n (x) := 2<M<?~ n , q n+1 ab; qa; q, q n x) - i</>i(0; qa; q, x) 

00 ( -\\k n k(k-i)/2 „k / k \ 

= E / w , ( - 1 + IK 1 " C 1 - 9" +J o6) J • 

^ (9a;?) fc (?;Q) fc V /J 7 

Since 

fc fc 
JJ(1_ Xj )>1_^ Xj if < X j < 1, j = 1, . . . ,k, 

we have 

k k k 

-1 + - Q n - J+1 ) (1 - q n+J ab) <J2<l n ~ J+1 + J2<l n+Jab 

3 = 1 3 = l J = 1 

g »-fc+i(l + a6) 

- 

if /c < n, while the left-hand side is equal to 1 otherwise. Thus, for |x| < M: 
, , N1 q n+1 (l + ab) ~ gfc(fc-i)/2 ? -fc M fc ~ q Hk-i)/2 M k 

Rn(x) <- ^ ^> ^ — ^ — + > -? — -, 

1-9 ^ {qa;q)k{q;q)k k ^ +1 {qa; q) k [q; q)k 

which tends to as n — > oo. □ 
Proposition A. 2. For a > — 1 we have 

lim 1 <j> 1 (0;q a+1 ;q,(l-q) 2 z)= F 1 (-,a + l;-z), 

uniformly for z in compact subsets of C. 
Proof. 

00 ( -i\fc _fc(fc-l)/2 /i _ „^2fc fc 

and the summand in the sum at the right-hand side can be majorized by 

(A f fl -«/*izn* TT (q(q+jV2 ~ ? 1+(Q+J)/2 ) (^ 2 - 9 1+ ^ 2 ) 
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Now, by [11, Lemma A.l] (read — 1 < [x — A instead of < [i — A in the formulation 
of that lemma) , we see that 



q (a+j)/2 _ q l+(a+j)/2 



1 _ ql+a+j 

increases to (1 + a + j) -1 as g t 1 if ot + j > 0. So, the expression in (A.l), for 
1/2 < g < 1, is dominated by 

(2 a / 2 \z\) k 



(a + l) k k\ 



if a > 



and by 



I Z I ^ 

const. r — — if — 1 < a < 0. 

(a + l) fc fc! 



So the proposition follows by dominated convergence. □ 

Appendix B. Proof of Propostion 4.2 

First observe that u(z) := 2^1 (a, b; c; g, z) (c not a non-positive integer power of 
q) is analytic in z on C outside the set {q~ l \ I = 0, 1, 2, . . . }, on which points the 
function has possible poles. This follows immediately by rewriting the second order 
(/-difference equation for -u, given in [6, Exercise 1.13], as 

(B.l) u(z) = (z- l)" 1 ^ + q b )z - q - 1 - l)u(qz) + (q ' 1 - q a+b z)u(q 2 z)) . 

Indeed, u(z) is analytic in z for \z\ < 1. Then (B.l) shows that u(z) can be extended 
to an analytic function in z for \z\ < q~ x with possible pole in 1, and then by iter- 
ation to an analytic function for \z\ < q~ l with possible poles in 1, q -1 , . . . , q~ l+1 . 
Now use the identity (see [6, (4.3.2)]) 

(B.2) 2<Pi (o, 0; c; z) = — 20i (a, ag/c; ag/6; g, cq/{abz)) 

yc, o/CL, z, q/ z, q)oo 

(a,c/b,bz,q/(bz);q)oo ,/,,,,, /f h \\ 

+ (ca/b^q/z-^ 2<M6 ' bQ/C; bQ/a; ^' Cq/{ahz)) ' 

where a/b and 2; are not an integer power of g, where c and cq/(abz) are not a 
non-positive integer power of g, and where a, 6, c, 2 7^ 0. Now put z := q n+1 /a 
(n = 0,1,2,...) in (B.2). Then the first term on the right-hand side of (B.2) 
vanishes and after some rewriting of the coefficient in the second term we obtain 

(B.3) 2 Ma,b;c;q,q n+1 /a) = W^Z^k 6 -» ^(qb/c, b; qb/a; q, q~ n c/b), 

where initially a/b and a are not an integer power of g, c and q~ n c/b are not a 
non-positive integer power of g, and a, 6, c 7^ 0. By continuity we see that (B.3) 
more generally holds if c, qb/a, q/a, q~ n c/b are not a non-positive integer power of 
q and a, 6, c 7^ 0. The case of (B.3) with n a non-positive integer then follows by 
redefining the parameters in (B.3). Finally we can state: 
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Lemma B.l. Formula (B.3) is valid for integer n ifc, qb/a, q/a, q n+1 /a, q n c/b, c/b 
are not a non-positive integer power of q and a,b,c^ 0. 

Now consider (2.5). Initially this is proved (see [6, (1.4.1)]) for \z\, \b\ < 1 and 
c, az not a non-positive integer power of q. By analytic continuation (2.5) is more 
generally valid for z, 6, c, az not a non-positive integer power of q. 

Let us turn to the proof of Proposition 4.2. Proceed with (4.5) as written just 
before Proposition 4.2 and take into account the constraints for (2.5) as just men- 
tioned. Then we initially obtain, under the assumptions that \sxy\ < 1 and x 2 ,y 2 
are not non-positive integer powers of q, that (4.7) holds with the first expression 
taken on the right-hand side if n > m and q n ~ m s~ 1 xy~ 1 not a non-positive integer 
power of q, and that it holds with the second expression taken on the right-hand 
side if m > n and q m ~ n s~ 1 yx~ 1 not a non-positive integer power of q. However, 
by Lemma B.l, the two expressions on the right-hand side of (4.7) are equal to each 
other for all integer values of m, n if s~ 1 xy~ 1 , q n ~ m s~ 1 xy~ 1 , s~ 1 yx~ 1 , q rn ~ n s~ 1 yx~ 1 
are not a non-positive integer power of q. Hence, under these last constraints the 
left-hand side of (4.7) is equal to both expressions on the right-hand side for all 
integer values of m, n. Finally, if s~ 1 xy~ 1 , q n ~ m s~ 1 xy~ 1 are not a non-positive 
integer power of q but s~ 1 yx~ 1 or q rn ~ n s~ 1 yx~ 1 is a non-positive integer power of 
g, then by continuity the left-hand side of (4.7) is still equal to the first expression 
on the right-hand side. Similarly, if s~ 1 yx~ 1 , q m ~ n s~ 1 yx~ 1 are not a non-positive 
integer power of q but s~ 1 xy~ 1 or q n ~ m s~ 1 xy~ 1 is a non-positive integer power 
of q, then the left-hand side of (4.7) is still equal to the second expression on the 
right-hand side. 

References 

1. G. E. Andrews, q-Series: their development and application in analysis, number theory, 
combinatorics, physics, and computer algebra, Regional Conf. Ser. in Math., no. 66, Amer. 
Math. Soc, Providence, R.I., 1986. 

2. G. E. Andrews and R. Askey, Enumeration of partitions: The role of Eulerian series and 
q-orthogonal polynomials, Higher combinatorics, M. Aigner, ed., Reidel, 1977, pp. 3-26. 

3. H. Exton, A basic analogue of the Bess el- Clifford equation, Jnanabha 8 (1978), 49-56. 

4. , q-Hypergeometric functions and applications, Ellis Horwood, 1983. 

5. Ph. Feinsilver, Elements of q-harmonic analysis, J. Math. Anal. Appl. 141 (1989), 509-526. 

6. G. Gasper and M. Rahman, Basic hypergeometric series, Cambridge University Press, 1990; 
second corrected and extended edition 2004. 

7. W. Hahn, Beitrage zur Theorie der Heineschen Reihen, Math. Nachr. 2 (1949), 340-379. 

8. , Die mechanische Deutung einer geometrischen Differenzgleichung, Z. Angew. Math. 

Mech. 33 (1953), 270-272. 

9. E. Heine, Handbuch der Kugelfunctionen, Theorie und Anwendungen. Erster Band: Theorie 
der Kugelfunctionen und der verwandten Functionen, Reiner, Berlin, 1878. 

10. M. E. H. Ismail, The zeros of basic Bessel functions, the functions Jv+ax(x)> and associated 
orthogonal polynomials, J. Math. Anal. Appl. 86 (1982), 1-19. 

11. T. H. Koornwinder, J acobi functions as limit cases of q-ultraspherical polynomials, J. Math. 
Anal. Appl. 148 (1990), 44-54.. 

12. M. Rahman, An addition theorem and some product formulas for q-Bessel functions, Canad. 
J. Math. 40 (1988), 1203-1221. 



20 



TOM H. KOORNWINDER AND RENE F. SWARTTOUW 



13. , A note on the orthogonality of Jackson's q-Bessel functions, Canad. Math. Bull. 32 

(1989), 369-376. 

14. L. L. Vaksman and L. I. Korogodskh, An algebra of bounded functions on the quantum group 
of the motions of the plane and q-analogues of Bessel functions, Soviet Math. Dokl. 39 (1989), 
173-177. 

15. G. N. Watson, Theory of Bessel functions, Cambridge Univ. Press, 1944. 

16. M. E. H. Ismail, D. R. Masson and S. K. Suslov, The q-Bessel function on a q-quadratic grid, 
Algebraic methods and ^-special functions, CRM Proc. Lecture Notes, 22, Amer. Math. Soc, 
1999, pp. 183-200. 

17. F. H. Jackson, The application of basic numbers to Bessel's and Legendre's fuuctions, Proc. 
London Math. Soc. (2) 2 (1904), 192-220. 

CWI, Amsterdam 

Current address: University of Amsterdam, Korteweg-de Vries Institute, P.O.Box 94248, 
1090 GE Amsterdam, The Netherlands 

E-mail address: T.H.Koornwinder@uva.nl 

Faculty of Technical Mathematics and Informatics, Delft University of Tech- 
nology 

Current address: VU University Amsterdam, Faculty of Exact Sciences, Department of Math- 
ematics, De Boelelaan 1081, 1081 HV Amsterdam, The Netherlands 
E-mail address: rene@few.vu.nl 



